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Abstract: We re-examine spontaneous CP violation at the tree level in the context 
of the next-to-minimal supersymmetric standard model (NMSSM) with two Higgs 
doublets and a gauge singlet field. We analyse the most general Higgs potential with- 
out a discrete Z% symmetry, and derive an upper bound on the mass of the lightest 
neutral Higgs boson consistent with present experimental data. We investigate, in 
particular, its dependence on the admixture and CP-violating phase of the gauge sin- 
glet field, as well as on tan (3. To assess the viability of the spontaneous CP violation 
scenario, we estimate ex by applying the mass insertion approximation. We find that 
a non-trivial flavour structure in the soft-breaking A terms is required to account for 
the observed CP violation in the neutral kaon sector. Furthermore, combining the 
minimisation conditions for spontaneous CP violation with the constraints coming 
from K°-K° mixing, we find that the upper bound on the lightest Higgs-boson mass 
becomes stronger. We also point out that the electric dipole moments of electron 
and neutron are a serious challenge for SUSY models with spontaneous CP violation. 



Keywords: Spontaneous Symmetry Breaking, Supersymmetric Models, Higgs 
Physics, CP Violation. 
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1. Introduction 

The origin of CP violation remains a fundamental open question in particle physics. 
In the standard model (SM), CP is explicitly broken at the Lagrangian level through 
complex Yukawa couplings which lead to CP violation in charged weak interactions 
via the Cabibbo-Kobayashi-Maskawa (CKM) matrix, V C km [!]■ Although the CKM 
mechanism can accommodate the experimental value of 6k (and, in principle, also 
s'/sk) in the neutral kaon sector, it is not clear whether it is the dominant contri- 
bution or the only one. An important motivation to consider new sources of CP 
violation stems from the fact that within the SM the amount of CP violation may 
not be sufficient to generate the observed baryon asymmetry in the universe [2]. 

An alternative scenario for the breaking of CP is to assume that it is a symme- 
try of the Lagrangian which is only spontaneously broken by the vacuum [3]. The 
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purpose of the present paper is to study spontaneous breaking of CP at the tree level 
within the context of supersymmetry (SUSY). Although the minimal supersymmet- 
ric standard model (MSSM) has two Higgs doublets, it is well known that due to 
SUSY constraints on the Higgs potential it is not possible to obtain spontaneous CP 
violation (SCPV) at the tree level. The reason for this is entirely analogous to the 
situation one encounters in non-supersymmetric two-Higgs-doublet models, where 
SCPV cannot be achieved if a Z 2 symmetry is imposed on the Lagrangian in order 
to guarantee natural flavour conservation in the Higgs sector [4]. The possibility 
that SCPV in the MSSM might be generated through radiative corrections has been 
explored in Ref . [5] , but this particular scenario has already been ruled out by exper- 
iment as it inevitably leads to the existence of a very light Higgs boson. Ultimately, 
this result is a consequence of the Georgi-Pais theorem [6] . It is therefore of interest 
to consider simple extensions of the MSSM such as a model with at least one gauge 
singlet field (JV) besides the two Higgs doublets (#1,2), the so-called next-to-minimal 
super symmetric standard model (NMSSM) [7,8], and to ask if one can achieve spon- 
taneous breaking of CP whilst generating the observed amount of ex and having 
Higgs-boson masses that are consistent with experimental data [9-13] 

Whether or not the observed CP violation in the kaon sector can arise solely from 
supersymmetry is a serious issue. This possibility has been examined in Refs. [14-16] 
by simply assuming Vckm to be real, and subsequently investigating whether SUSY 
sources of CP violation can account for the observed magnitude of €k and e'/ex- 
Here we would like to draw attention to the particularly attractive SUSY scenario 
with spontaneous CP symmetry breaking. The main point of this class of models 
is that the reality of the CKM matrix is automatic (see, e.g., Refs. [4,17]), and no 
longer an ad hoc assumption. In the specific scenario we shall be considering, CP 
violation is caused by the phases 0£> and 0tv associated with the vacuum expectation 
values of (H®) and (N) respectively. Being a singlet Higgs field, N does not couple 
to quarks, and although does couple to the up-type quarks, the phase appearing 
in the quark mass matrix can be rotated away by redefinition of the right-handed 
quark fields. Consequently, this phase does not appear in the CKM matrix since 
the lU-boson interactions are purely left-handed, but it does enter in other SUSY 
interactions. 

The plan of this paper is as follows. In Sec. 2, we analyse the conditions required 
for the Higgs potential to have a CP-violating global minimum. We show that the 
upper bound on the mass of the lightest Higgs boson imposes constraints on the 
CP-violating phases, ipo and <pN- Section 3 is devoted to a brief review of the real 
CKM matrix in supersymmetric extensions of the SM with two Higgs doublets and 
an arbitrary number of gauge singlet fields. In Sec. 4, we introduce the mass and 
mixing matrices, as well as the couplings of the Higgsino, lU-ino, and singlet field in 
the NMSSM with spontaneous CP violation. The calculation of the relevant SUSY 
contributions to e K in the mass insertion approximation is presented in Sec. 5. We 
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pay particular attention to the case of low tan (3 and examine the question of whether 
the choice of parameters of the Higgs potential that is required to obtain spontaneous 
CP violation is consistent with the experimentally observed indirect CP violation in 
the kaon sector. The implications of for the upper bound on the lightest Higgs- 
boson mass and a new flavour structure, besides the one of the CKM matrix, are 
discussed. We also comment on the electric dipole moments of electron and neutron. 
Finally, we summarize the key features of our analysis in Sec. 6. 

2. The Higgs potential 

2.1 The superpotential 

We are concerned with the next-to-minimal supersymmetric standard model dis- 
cussed by Davies et al. [18]. The most general form of the superpotential that we 
analyse here is 1 

W = Wfermion + W^ffiggs, (2.1) 

with 

W^fermion = S ab (h^UjH^ + h%Q\D jH? + ^RjH^ , (2.2) 

W^Higgs = -\Ne ab H«H b 2 — ^ iV 3 — rN — fie ab H?H b 2 , (2.3) 

where i,j = 1,2,3 denote generation indices, a, b = 1,2 are SU(2) indices, and e is 
a completely antisymmetric 2x2 matrix with e 12 — l. In the above expression, H± 
and H 2 denote the Higgs doublets of the minimal supersymmetric standard model 
and N is a singlet field. The matrices hu,ho, and He give rise to the usual Yukawa 
interactions which generate the masses of quarks and leptons. Since we are dealing 
here with spontaneous breaking of CP rather than explicit CP violation in the Higgs 
sector, these matrices are real. 

2.2 Soft-breaking terms 

In addition to the superpotential given by Eq. (2.1), we have to specify the SUSY 
soft-breaking terms. These are given by 

£ S B = 4 C B mi ° n + 45^, (2-4) 

where 

+ M 2 Ri .R-R) - \ (M s \ s \ s + MA A + M'X'X' + H.c.) 

+ e ab ^h^U 3 H b 2 + AthtQlD^ + A%h%L\R 3 H«) , (2.5) 

1 This form of the superpotential has also recently been studied in Rcf. [19], within the context 
of electroweak baryogenesis. 
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and 



s = m \ x H?Hl + m 2 H2 H^E a 2 + m 2 N N*N 
- ^Bfxe ab H^H b 2 + A x Ne ab H«H b 2 + ^N 3 + A r N + H.c.) . (2.6) 

2.3 The scalar potential and spontaneous CP violation 

Following Ref. [18], we do not require the superpotential to be invariant under a 
discrete Z 3 symmetry which would imply /i — r — (thereby solving the so-called 
'yu problem' of the MSSM). 2 Further, we do not relate the soft SUSY-breaking pa- 
rameters to some common unification scale but rather take them as arbitrary at the 
electroweak scale. 

As it was noted by one of the authors (J.C.R.) a long time ago [12], the NMSSM 
with a Z3 symmetry has no spontaneous CP violation (no-go theorem). The inclu- 
sion of the Z 3 -breaking terms in Eq. (2.1), on the other hand, evades that no-go 
theorem, as was shown by Pomarol [13]. Throughout we shall assume that the tree- 
level potential is CP conserving and take all parameters real, but allow complex 
vacuum expectation values (VEV's) for the neutral Higgs fields which emerge after 
spontaneous symmetry breaking: 

From the superpotential and soft supersymmetry breaking terms, Eqs. (2.3) and 
(2.6), we derive the following CP-invariant neutral scalar potential: 

V = V + + X 2 vl) + i[(A 3 + X 4 )vfv 2 2 + (X 5 vl + \vl)vl\ 

1 11 

+ -X 7 viv 2 vl cos(6>i + 2 - 20 3 ) + -A 8 v 3 + -j=X/j,(v 2 + v 2 )v 3 cos(6> 3 ) 

+ \{rn\vl + m\v\ + rri%v\) -=m±vxv 2 v 3 cos(#i + 6 2 + 3 ) 

2 V 2 

22 



-.m 5 v 3 cos(3# 3 ) + m 6 viv 2 cos(#i + 6 2 ) + m 7 v 3 cos(26> 3 ) 



3V2 

+ V2mlv 3 cos(9 3 ), (2.8) 



where Vq = r 2 and 



Xi = \2 = \(g 2 + g' 2 ), \ 3 = -\(g 2 + g ,2 -4\ 2 ), A 4 = 0, 

A 5 = A 6 = A 2 , A 7 = Xk, X 8 = k 2 , (2.9) 

m i,2 = m H lt2 + /A m l = m /V! m 4 = A x , m 5 = A k , 

ml = —Bji + Ar, m 2 = kr, ml = —A r . (2-10) 



2 We remind the reader that a spontaneously broken discrete symmetry may lead to cosmological 
domain- wall problems [20]. We do not pursue this subject here. 
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In the numerical analysis reported here, we take the A;'s (I — 1, 2, . . . , 8) as given by 
Eq. (2.9), i.e. in terms of g, g', A, k, but consider the m^s as arbitrary parameters. 
Referring to Eq. (2.8), we see that the doublet phases enter only in the combination 
Q\ + 02, and thus we may define 

< p D = 6l + e 2j <p N = e z (2.ii) 

for the doublet and singlet phase respectively. We shall henceforth set 9\ = in 
Eq. (2.11), without loss of generality. 

2.4 Mass spectrum 

To obtain the physical masses of the NMSSM with SCPV, the following procedure has 
been adopted. We randomly choose sets of parameters and perform the minimisation 
numerically. Next, we check that the local minimum is a true one not only by 
verifying the minimisation conditions, but also by ensuring that all the physical 
Higgs bosons have positive squared masses. In calculating the neutral Higgs-boson 
mass matrix, radiative corrections due to top-stop loops have been included. Finally, 
we require that the chargino and squark masses are above the current experimental 
lower limits [21]. 

We find that an acceptable mass spectrum can be easily obtained, with the exact 
values depending on the set of parameters we choose. Rather than presenting tables 
or plots for different sets of parameters, we will concentrate on two issues: (a) the 
influence of the CP-violating phases on the lightest Higgs-boson mass; and (b) the 
maximum theoretically attainable value of its mass. As far as the former is concerned, 
we confirm the results of an analysis of the Higgs potential performed in Ref. [18], 
where it was pointed out that the large-phase solution is favoured. This feature can 
be clearly observed in Fig. 1, where we display the Higgs-boson mass, m#o, as a 
function of the phase of the gauge singlet field. It is important to note that large CP 
phases may be in conflict with existing limits on the electric dipole moments (EDM's) 
of electron and neutron. However, these constraints can be evaded, for instance, if 
substantial cancellations among different SUSY contributions do occur [22,23] or if 
non-universal soft-breaking terms are present [24]. (We will return to this point in a 
subsequent section.) 

As for the maximal possible value of the Higgs-boson mass, the result can differ 
from that of the MSSM for the case of large values of the coupling constant A. In 
Fig. 2, we show the tree level and one-loop corrected Higgs-boson mass as a function 
of A. It is apparent that only for large values of the coupling A the situation is 
different from that of the MSSM. For low values of A, corrections to the tree level 
Higgs-boson mass are significant and depend mainly on the SUSY scale that we take 
for the squarks. This is illustrated in Fig. 3 where we plot the upper limit on m#o 
versus tan/3 = v 2 /v 1 [cf. Eq. (2.7)] for different values of M SUSY . Before leaving the 
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Figure 1: Maximum value of the lightest Higgs-boson mass (in GeV) as a function of the 
CP-violating phase 4>N (in radians) of the Higgs singlet field. 
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Figure 2: Maximum value of the lightest Higgs-boson mass (in GeV) vs A at the tree level 
and after including radiative corrections (at one-loop level) for Mstjsy = ITeV. 



subject of the Higgs potential, it should be emphasized that the Higgs boson mass 
limits obtained at LEP for the standard model and the MSSM (see, e.g., Ref. [25]) 
do not necessarily apply to the NMSSM because the lightest neutral Higgs boson 
may have a reduced coupling to the Z° due to some singlet admixture [7,26]. In this 
case, the Higgs boson might not be detectable. 

3. A real CKM matrix 



We outline briefly the scenario of a real CKM matrix in the framework of supersym- 
metry with SCPV. In this case, CP invariance is imposed on the Lagrangian, and 
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Figure 3: Upper bound on the lightest Higgs-boson mass (in GeV) as a function of tan/? 
for various values of Mstjsy- 

consequently all couplings are real. To address the question of whether the CKM 
mechanism occurs in the NMSSM model, let us consider the quark Yukawa couplings. 
As mentioned in the introduction, the gauge singlet field does not couple to quarks. 
Thus, the Lagrangian in terms of the weak eigenstates may be written as (omitting 
generation indices) 

C Y = -huH^ULUR - h D H°*d L d R + H.c., (3.1) 

where hu,D are arbitrary real matrices in flavour space. After spontaneous symmetry 
breaking, the up- and down-type quarks acquire masses, namely 

mu = hu^=e- i<t,D , m D = h D ^=. (3.2) 

Note that the overall phase <pD can be rotated away by means of a phase trans- 
formation on ur, i.e. ur — * u' R = e~ t ^ >D UR. Since the VT-boson interactions are 
purely left-handed, this phase does not show up in the charged weak interactions. 
Consequently, the CKM matrix is real and CP violation arises solely from the rel- 
ative phases that appear in the VEV's of the neutral Higgs fields, Eq. (2.7). An 
alternative way to derive the above result is to compute the weak basis invariant 
T = Tr[Hu, Hp} 3 , where Hy = mum\j and Hp = m^m^ [27]. It follows from 
Eq. (3.2) that T vanishes, and hence there is no CP violation through the CKM 
mechanism in our model. 

4. Content of the NMSSM 

As we have argued in the previous section, all supersymmetry CP-violating phases 
are equal to zero and the CKM matrix is real. Even so, the phases associated with the 
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VEV's, <pD and </>tv, appear in the scalar quark, gaugino and Higgsino mass matrices, 
as well as in some of the vertices. We first consider the squark mass matrices. 

4.1 Squark mass matrices 

The 6x6 squark mass-squared matrices in the (qL, Qr) basis are given by 

M- = ( M ^ LL M ^ LR \ a = U D (4 1) 

\ 1RL QRR / 

made up of 3 x 3 submatrices. It proves convenient to work in the so-called 'super- 
CKM' basis, in which the quark mass matrices are diagonal, and squarks as well 
as quarks are rotated simultaneously (see Ref. [28] for details). In this basis, the 
aforementioned submatrices are of the form 

MI ll = (mj iag ) 2 + VEMlV^ + \m\ cos 2/3(3 - 4 sin 2 6 w )t, 

M kn = K 1 ^) 2 + V n M k V R + l M z cos ^ sin2 
MI lr = Mil = VEY*V^^ - cot Pe^m^, 



M% ll = (m^) 2 + V?M%V* - l -Ml cos 2/3(3 - 2 sin 2 6 w )t, 
M kn = «*? + V n M k V n' ~ cos2/?sin 2 ^l, 



M ku = M L = V L DY n V R^ - ^tan/fe^m** (4.2) 



6w being the Weinberg angle, 1L denotes a 3 x 3 unit matrix, and the \^'s diago- 
nalize the up- and down-type quark mass matrices in Eq. (3.2). [Note that the V 
matrices are orthogonal, and therefore Vckm = Vl^Vl ) T -] m what follows, we take 
a particular case of the 'super-CKM' basis where Vf = 1 and = V'ckm, so that 



? s ) 2 = VcKM^m^V^KM. (4.3) 



Further, we have introduced the shorthand notation 

Yg j = A^h^, q = U,D (no sum over (4.4) 
fjLeS = pL + \^ e i*N (4_ 5 ) 

[Notice that Eqs. (4.2)-(4.5) show explicitly the dependence on the CP phases.] As 
we shall see in the next section, a non-universal flavour structure in the A terms, i.e. 
A l i 7^ constant, is indispensable for having sizable supersymmetry contributions to 
CP violation in the kaon sector. 3 At this stage, we therefore make no assumptions 
regarding the nature of the several matrices and couplings involved. 

3 For related work, see Refs. [14, 15]. 
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The squark mass-squared matrix, M|, can then be diagonalized by a 6 x 6 unitary 
matrix such that 

{Mf as ) 2 = R^MlRt. (4.6) 

4.2 Chargino mass matrix 

We now turn our attention to the four-component Dirac fermions, Xi,2-> which arise 
from the mixing of the H^-inos and the charged Higgsinos. Within SUSY, the 
chargino mass terms in the Lagrangian are 



2 v M x o ; \^ 

where i\) +T = (—iX^.H^)-, ^>~ T = (—i\~,Hi). The chargino mass matrix in the 
NMSSM with spontaneous CP violation reads 

M = ( Ms V2M w sm(3e-^\ 

x \V2M w cosp fi cS J 1 ' ; 

with the gaugino mass M 2 being real and /j, e g given in Eq. (4.5). As usual, can 
be diagonalized by means of a biunitary transformation, i.e. 

M diag = u* M _y^ (4_ 9 ) 

where M? iag is positive and diagonal, and U, V are unitary matrices. 

For convenience, we perform our calculation in the weak (rather than the phys- 
ical) basis, i.e. W~ = (—i\~,i\ + ) T , H~ = (H[ , H 2 ) T , and their relation can be 
summarized as follows: 4 

PlW- = Pl(UZ 1 xI + UZi&), PrW- = P r (V 11 xT + V 21 &), (4.10a) 

P l H- = Pl(U: 2 xI + U; 2 x;), PrH- = Pr(V 12 xI + V 22 x;), (4.10b) 

with Pl,r = (l=F75)/2. In order to specify the relevant couplings within the NMSSM, 
we recall that [29] 

-£^,H = g (HfwP L H + H°*HP L W) + XNHP L H + H.c. (4.11) 

Then, substituting the VEV's from Eq. (2.7), and adding VY-ino and Higgsino 'mass' 
terms, we arrive at 



-C mt = m w WW + nifjHH + -^(v^W R H L + v 2 e^W L H R + H.c), 

a/2 



(4.12) 



4 We define W and H as particles, contrary to Ref. [29]. 
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where 

m w = M 2 , rrijj = |/x cfT |, ip = arg (/x efT ) , (4.13) 

and ii e fi as in Eq. (4.5). In deriving Eq. (4.12), we have adjusted our definition of the 
phase of the left-handed Higgsino field such that Hl — > e~ %Lp H.L. As for gluino and 
neutralino interactions, we will argue in the next section that to good approximation 
the chargino box diagram may be regarded as the dominant contribution to €k (see 
also Refs. [11,13]). 



5. Implications of indirect CP violation for the NMSSM 

In this section, we wish to explore the consequences of SCPV for the upper bound of 
the lightest Higgs-boson mass derived in Sec. 2 by taking into account CP violation in 
mixing. To accomplish this, we will compute the box-diagram contributions 
to tK by applying the mass insertion approximation [30]. That is, adopting the 
notation 

Mf = (M°) 2 + Mf, (M°) 2 = diag(a?,...,a°), (5.1) 
and recalling Eq. (4.6), we may write to first order in M~ 2 [31] 

(R\)mkf(a k )(Rg) kn = Smnf(a° n ) + (Mf) mn F(a° m ,a° n ), (5.2) 
/ being an arbitrary function, k — 1, ... ,6, and 

F(l , y) = /(*>-/(»> , (5.3) 
x-y 

5.1 Effective Hamiltonian 

Let us start with the effective Hamiltonian governing AS* = 2 transitions, which can 
be written in the form [32]: 

H cS = J2c l O i , (5.4) 

i 

where the q's are the Wilson coefficients describing the short-distance interactions 
while Oi denote local operators which may be found, e.g., in Refs. [33,34]. The 
off-diagonal element of the neutral kaon mass matrix is related to H ef[ through 

Ma = 

2m K 

and its imaginary part gives rise to (assuming e' €k) 

e ff/4 ImM 12 

€K * ^"A%' (5 - 6) 
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Figure 4: The dominant box-diagram contribution to the off-diagonal element ImA^i2 in 
the neutral kaon mass matrix within the framework of the NMSSM. 

with the experimental values of Am K = (3.489 ± 0.008) x 10~ 12 MeV 5 and \e K \ = 
(2.271 ±0.017) x lO" 3 [21]. 

The Wilson coefficients q in the presence of SUSY contributions can be decom- 
posed as follows: 

Ci = cf + cf +cf + c? + cf. (5.7) 

Regarding the various contributions in Eq. (5.7), we confine ourselves to a relatively 
few remarks. First, there are no VF-boson and charged Higgs-boson contributions 
to the imaginary part of M.12 since the CKM matrix is real. Second, in subsequent 
calculations we will choose a basis where the down-type Yukawa matrices are diag- 
onal. As a consequence, the flavour off-diagonal gluino contributions are either zero 
or real in the approximation of retaining only a single mass insertion in an internal 
squark line. Third, box diagrams with neutralinos are expected to compete with the 
chargino contributions only for large values of tan/5 (i.e. tan/3 ~ 50), and hence are 
less important for the present scenario with low tan (3. In the remainder of this sec- 
tion, we will therefore focus on the chargino box diagram with scalar up-type quarks 
(see Fig. 4) as it gives by far the dominant contribution in our model. 

Turning to the operators Oi in Eq. (5.4), the AS* = 2 transition within the 
model under study is largely governed by the V-A four-fermion operator 0\ = 
d^Phsd^^PLS. As a matter of fact, the new-physics contribution to the remain- 
ing operators (i.e. scalar and pseudoscalar operators) are suppressed relative to the 
chargino contribution either through small quark masses, m St d, or small Yukawa cou- 
plings. Thus, it is fairly reasonable to consider only the non-standard contribution to 
the Wilson coefficient C\. The relevant hadronic matrix element in the vacuum sat- 
uration approximation is given by (K°\Oi\K°) = (2/3)/|-m|-, with f K = 160 MeV. 

5.2 Chargino box diagram 

In the W-H basis, and making use of the mass insertion approximation, the chargino 

5 Note that for Amjj, and following standard procedure, we have used the experimental value. 
No useful bound on the present model can be derived from Am^| oxp , since it receives important 
long-distance contributions. 
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Figure 5: The main contributions to €k in the mass insertion approximation with IF-ino 
and Higgsino exchange. 



contributions to the short- distance coefficient c\ can be written as a sum of 16 in- 
dividual terms, each corresponding to a particular box diagram. Introducing the 
notation {ij} = L, R) for the various mass insertions, these contributions can be 
classified in terms of the nature of the mass insertion in the internal squark line and 
the number of W-H exchanges in the sfermion line. That is, the relevant amplitude 
for ci can be symbolically written as 

~ (X W ^)°[{LL}{LL} + {RR}{RR} + {LR}{RL}] 
+ {Xw-hY[{LL}{LR} + {LR}{RR}] 

+ (Xw-h) 2 [{LL}{RR} + {LR}{LR}\. (5.8) 

Only one diagram belongs to the {LL}{LL} and {RR}{RR} mass insertions, whereas 
we find two (four) diagrams associated with the {LR}{RL}, {LL}{RR}, {LR}{LR} 
({LL}{LR}, {LR}{RR}) mass insertions. It turns out that the dominant contribu- 
tion to Ci corresponds to the class of diagrams with (X^^^-fLLj-fLi?}, depicted 
in Fig. 5, in agreement with the findings of Refs. [11, 13]. The calculation of the 
dominant contributions to Im7Wi2 for the case of non-degenerate left-handed up- 
squarks leads to a cumbersome expression. Hence, for the sake of simplicity, we 
perform the calculation in the limit of degenerate left-handed up-type squarks, and 
find that the expression for Im.Mi2 becomes much simpler, and its physical inter- 
pretation more transparent. (We have numerically verified that small values of the 
off-diagonal elements of M~ have little effect on the ul mass splitting.) In that 
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limit, i.e. rriu L = (rriq), where (m,-) represents the average squark mass of the first 
two generations, we obtain 

lmMi2 = 2Glflm K ml + ^ ^ 

OH \ m q) 

x {AA uS m(<p x -(f) D ) (M?) i2 /(r^,r^,r SL ,r 4fl )}, (5.9) 
with rrifj ^ and (M?) as defined in Eqs. (4.13) and (4.2) respectively. Here we have 

' Q 12 

used unitarity of the CKM matrix and kept only the leading top-quark contribution. 
In the above formula, r\ = mf/ (rriq) 2 , AAu = A]f — A^, and 

ip x = arg {e i4>D my V + cot (3m jj) ■ (5.10) 

The function / in Eq. (5.9) can be reduced successively to elementary functions for 
appropriate limits of its arguments (see Refs. [23,31]). For ease of writing, we give 
the result for the special case 



I(x,x, 1,1) = i 



10 + 19x + a; 2 (l + 6x + 3x 2 ) 
3(1 -x) 5 + (1 -x) e nX 



(5.11) 



which corresponds to the scenario with = and degenerate squarks. The 
appearance of Ay (i = 1,2) terms in Eq. (5.9) is closely related to having the 
Higgsino coupling to down and strange quarks in distinct diagrams, and stems from 
the trilinear soft-breaking terms in Y$ = A^hy [Eq. (4.4)]. From inspection of 
Eq. (5.9), it is straightforward to conclude that in order to get Im different from 
zero we need to move aside from a theory of universal A v terms (i.e. AAu 7^ 0) - 
in other words, it is not possible to saturate the observed CP violation in the K- 
meson system in the context of SUSY with a real CKM matrix and universal Ajj 
terms [14, 15]. 6 Let us note parenthetically that the need for a special A§ (i = 1, 2) 
texture as a key element to get the experimentally measured €k has also been pointed 
out in Ref. [14], although in a different scenario. 

Since the sign of AAu is not fixed, we do not include the constraint Re ex > 
in the analysis which follows. 

5.3 Numerical results and discussion 

For our numerical calculations, we have used the nominal values 

(M| )J (rriq) 2 = 0.08, V td = 0.0066, V ts = -0.04, m t = 175 GeV. (5.12) 

Notice that the value of V t d differs from the one determined in the context of the SM 
with CP violation by fitting to experimental data; this is due to the fact that we are 

6 It has recently been pointed out in Ref. [35] that there might be a connection between a non- 
trivial flavour structure and sizable SUSY CP violation, thereby avoiding domain-wall problems. 
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Table 1: Numerical values of |e#-| in the low tan/3 region for certain sets of model pa- 
rameters that satisfy the minimisation condition of the Higgs potential. We have chosen 
AAjj = 500 GeV for the non-universal A terms, as described in the text. 



\fjA 


WD 






( in ~ ) 
\" b ql 


" l tR 


tan 3 


A 


v i 


(10 ) 


(racl ) 


(racl ) 


(GeV) 


(GeV) 


(GeV) 






(GeV) 


3.24 


4.71 


1.57 


99 


252 


235 


6.7 


-0.03 


327 


3.03 


0.89 


1.75 


97 


261 


168 


6.6 


+0.33 


387 


2.75 


4.71 


4.71 


99 


232 


201 


9.2 


-0.02 


221 


2.42 


1.96 


4.08 


94 


299 


174 


5.1 


-0.06 


352 


2.10 


4.67 


4.75 


98 


279 


220 


7.8 


+0.01 


142 


2.02 


4.68 


4.71 


92 


250 


152 


7.4 


+0.02 


371 


2.01 


4.18 


4.73 


96 


280 


232 


4.6 


-0.01 


238 


1.31 


1.12 


4.72 


100 


273 


241 


9.6 


-0.01 


238 


1.29 


2.35 


4.70 


99 


258 


230 


6.1 


-0.13 


363 



dealing with a flat triangle and the corresponding orthogonality constraints. 7 The 
remaining parameters are evaluated numerically by minimising the Higgs potential 
(see Sec. 2). Our results for the absolute value of ex for various sets of SUSY 
parameters and low tan/5 are reported in Table 1. Recall that the numerical values 
we present for |ex| have been derived for degenerate left-handed squarks. Taking into 
account corrections due to non-degeneracy will lead to a very small enhancement, 
typically a few percent, for Im.Mi2- Hence, the values for \e K \ listed in Table 1 may 
be interpreted as a lower bound. As far as M 2 and /i are concerned, we obtain the 
ranges 90 GeV < M 2 < 160 GeV and 120 GeV < < 270 GeV. Turning to the 
non-universal flavour structure, it is obvious from Eq. (5.9) that there is a linear 
dependence of €k on the relative difference AAjj. In order to saturate the value of 
| | and to obey present experimental limits on the sparticle spectrum, one has to 
take AAu of order 500 GeV. A detailed discussion of scenarios where such values 
for Ay (i = 1,2) appear lies beyond the scope of this paper. We will just mention 
that values around the TeV scale do not significantly affect the mass spectrum of the 
theory, and that they can account for values of the left-right mass insertion, i.e. 

(fV \ _ V ? ( V L Y U V R )«3 ( r-,o\ 

V2 (rriq) 

as small as 0.02 and 0.08 for (5^ R ) 13 and (5^ R )23 respectively, which are consistent 

7 The value used for \Vtd\ is consistent with the experimental value of B^-B^ mixing, provided 
one assumes a new-physics contribution to Ams d of at least 20%. To be specific, we have chosen 
a positive value for p. As for a negative value of p, orthogonality of the CKM matrix would imply 
\Vtd\ =0.011. 
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with the bounds coming from the measurement of b — > S7 branching fraction. (Note 
that these <5's also affect the decay b — > sl + l~.) Indeed, according to Ref. [28], 
(5^)23 is constrained to be 1(^^)23! ^ 3({ m q) /500GeV) 2 , which is only useful if 
{rrig) < 300 GeV. Similarly, the chargino contribution to B®-B® mixing leads to the 
constraint (S^ R ) 13 < 0.1((m ? ~) /500 GeV) [31]. 8 

At this point, a few remarks are in order regarding previous studies contained in 
the literature. As we have already mentioned, the diagrams depicted in Fig. 5 give 
the leading contributions to Im.Mi2, in agreement with Refs. [11,13]. The final form 
of our result, however, differs from that of previous works in several ways. First, we 
have used the mass insertion approximation to compute the dominant contributions 
to ex, which enables us to use other largely model-independent constraints on the 
off-diagonal elements of the squark mass matrices. Second, we have considered a non- 
trivial structure for the LR squark mass-squared terms, which proved to be crucial 
in order to have indirect CP violation within the present model. More importantly, 
our parameter space is severely constrained since in the calculation of the various 
contributions to ex we have considered only those sets of parameters that correspond 
to true minima of the Higgs potential. 

We conclude this section with several comments on large CP phases and their 
implications for the EDM's of electron and neutron. 

5.4 Remarks on the electric dipole moments 

As seen from Table 1, the CP phase of the gauge singlet field, 4>n, is accompanied 
by small values of A so that the EDM constraints are less stringent [37]. Further, 
the SUSY contributions to the EDM's depend only on flavour-diagonal entries in 
the squark mass matrix, whereas the expression for €k involves flavour off-diagonal 
elements. However, it is evident from Eq. (4.2) that there is an overall CP-violating 
phase which is independent of the new flavour structure. Apart from the requirement 
that A l jj\ i= j <C Ayli-^j, a conceivable way to suppress the EDM's includes large 
cancellations between different SUSY contributions and a SUSY particle spectrum 
with masses of scalar quarks and gauginos in the TeV range. Given the analytic 
results for the contributions to the EDM's of electron and neutron mediated by 
photino and gluino [33], together with the sets of parameters displayed in Table 
1 and the present experimental results of d n < 6.3 x 10~ 26 ecm (90% C.L.) and 
d e = 1.8 x 10 _27 ecm [21], the photino and gluino masses are required to satisfy 
0.5 TeV <m;<2 TeV and 2 TeV < m~ q < 6 TeV. 9 

Such a hierarchy in the soft gaugino masses requires further discussion. We first 
note that masses of the superpartners of about 1 TeV may be in conflict with the 

8 Note that the requirement of vacuum stability also leads to upper bounds on (^ fl )i3 (i = 1, 2) 
[36]. The nominal values used in our numerical analysis are compatible with these upper limits. 

9 For simplicity, wc have assumed mj — nig. Note that the constraint on m g is relaxed if we 
allow for an order-of-magnitude variation in the value of the hadronic matrix elements. 
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cosmological relic density. In addition, such a scale for the 7 and the g masses is 
rather unnatural, since the masses of the squarks and lU-ino are typically of the order 
100-300 GeV in the model under consideration. [Similar conclusions can be drawn if 
we take into account the chargino and neutralino contributions (see, e.g., Ref. [23]).] 
Finally, we stress that the above-mentioned hierarchy for the scalar partners leads 
to an unacceptable scenario for the lightest supersymmetric particle (LSP). In this 
case, the LSP would be either charged (one of the lightest squarks or a chargino), 
or would have a non-zero lepton number. To summarize, it appears that unless 
the parameters of the model are fine-tuned, there will be large contributions to the 
EDM's of electron and neutron, in potential conflict with experiment. 10 

6. Conclusions 

The origin of CP violation, in particular, the question of whether CP is explicitly or 
spontaneously broken, is an important issue in particle physics. In this paper, we have 
studied spontaneous CP violation in the context of the NMSSM, demonstrating that 
it is possible to generate sufficient CP violation in order to account for the magnitude 
of e K . We have emphasized that the NMSSM with spontaneous CP violation is a 
natural framework to discuss the possibility of having CP violation that arises entirely 
from SUSY interactions. This is due to the fact that in the above-mentioned scenario 
the CKM matrix is automatically real. 

By performing a complete and systematic study of the Higgs potential in the 
NMSSM with spontaneous CP violation, we have shown that the minimisation of 
the most general Higgs potential leads to an acceptable mass spectrum which is 
accompanied by large CP- violating phases. We have argued that the lightest neutral 
Higgs boson can have a reduced coupling to the Z° due to the additional singlet field, 
and thus may escape detection at LEP. We have shown that it is possible to have a 
mass of the lightest neutral Higgs boson of about 125 GeV. However, such values for 
the Higgs-boson mass require a relatively large value for M S usy of the order 0(1 TeV). 
By contrast, a rather low SUSY scale with M S usy ~ 300 GeV (i.e. light squark and 
PU-ino masses) and a non-trivial flavour structure of the soft SUSY-breaking trilinear 
couplings A l {j (i = 1,2) are required in order to account for the observed CP violation 
in K°-K° mixing. As a consequence, the parameter space is severely constrained and 
the mass of the lightest Higgs boson is further diminished, and it turns out to be no 
greater than ~ 100 GeV for the case of low tan/9 (< 10). As far as large CP phases 
are concerned, we have argued that it may be difficult to reconcile the large-phase 
solution with the severe constraints on the EDM's of electron and neutron. Although 
we do not exclude the possibility of cancellations between different contributions, we 

10 As far as the B system is concerned, we merely remark that the model can accommodate a 
large CP asymmetry aj/^, i.e. a large value of sin 2/3 (see, e.g., Ref. [10]). 
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find this highly unlikely in view of the constrained parameter space. Therefore, the 
implications of the EDM bounds on the parameter space, as well as the implied LSP 
scenario, are a great challenge for SUSY models with spontaneous CP violation (at 
least within the minimal model we have considered here). 
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